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Abstract. IceCube is currently not only the largest neutrino telescope but also one of the world's most competitive instruments 
for studying cosmic rays in the PeV to EeV regime where the transition from galactic to extra-galactic sources should 
occur. Further augmenting this observatory with an array of radio sensors in the 10- 100 MHz regime will additionally permit 
observation of the geomagnetic radio emission from the air shower. Yielding complementary information on the shower 
development a triple-technology array consisting of radio sensors, the ground sampling stations of IceTop and the in-ice 
optical modules of IceCube, should significantly improve the understanding of cosmic rays, as well as enhance many aspects 
of the physics reach of the observatory. Here we present first results from two exploratory setups deployed at the South Pole. 
Noise measurements from data taken in two consecutive seasons show a very good agreement of the predicted and observed 
response of the antennas designed specifically for this purpose. The radio background is found to be highly dominated by 
galactic noise with a striking absence of anthropogenic radio emitters in the frequency band from 25-300 MHz. Motivated 
by the excellent suitability of the location, we present first performance studies of a proposed Radio Air-Shower Test Array 
(RASTA) using detailed MonteCarlo simulation and discuss the prospects for its installation. 
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I. INTRODUCTION 

Detection of air- showers through geomagnetic radio 
emission in the MHz regime offers various advan- 
tages compared to other detection technologies. In con- 
trast to ground-based shower sampling, the radio sig- 
nal represents the integral shower development and 
hence promises excellent energy resolution. In contrast 
to fluorescence measurements, high duty-cycles can be 
achieved with a low-cost technology. Further, the lateral 
distribution of the radio signal will provide an indepen- 
dent measure of the depth of the shower maximum. Con- 
sequently, radio extensions are now being developed in 
many existing air-shower experiments. 

The IceCube observatory at the South Pole employs 
a 1km 3 volume, instrumented with an array of photo- 
sensors as an ice-Cherenkov detector mainly targeted at 
the detection of astrophysical neutrinos. Installed at a 
depth between 1 450-2450 m, it is also sensitive to muons 
from air- showers above an energy at the surface of ~ 
300 GeV. This detector is complemented by a 1 km 2 sur- 
face array of ice-Cherenkov tanks that sample electrons 
and muons with a threshold of lOMeV from air-showers. 
Coincident detection of high-energy muons in ice and the 
shower particles on the surface provides a unique mech- 
anism to study cosmic ray composition from the knee to 
the ankle with a primary threshold of ~ 300 TeV, albeit 
with a limited aperture of 0.3 km 2 sr. 



Extending this observatory with a large array of radio 
antennas can significantly enhance the physics potential 
of this observatory [1]: adding a third detection technol- 
ogy with complementary shower observables will per- 
mit significant reduction of the systematic uncertainties 
in the composition measurement. Extending to several 
km 2 , it may serve as a veto for cosmic ray air- showers, 
and thus enhance the neutrino sensitivity in the southern 
hemisphere. Inverting the veto logic, it may also serve 
as an instrument for ultra-high energy gamma-rays by 
searching for muon-poor air- showers. 

In this paper, we show results from exploratory studies 
of the South Pole for the radio detection technique, as 
well as a first simulation of the performance that can be 
achieved with a dedicated test array. We also discuss the 
prospects for this approach. 

2. SITE EXPLORATION 

To establish the feasibility of radio air-shower detection 
at the South Pole beyond the proof of principle provided 
by the LOPES [2], Codalema[3] and AERA[4] experi- 
ments, two question applying to this particular location 
must be addressed: is it possible to build antennas that 
operate reliably at the environmental low temperatures; 
and does the ambient noise permit a sufficiently low trig- 
ger threshold. 
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FIGURE 1. A fat tube-dipole deployed at the South Pole. 
Trenches are backfilled after deployment. 

2.1. Antenna design 

The harsh environment at the South Pole, in particu- 
lar the low temperatures, pose significant challenges to 
antenna construction. On the other hand, the good radio 
transparency of the upper most layer of the glacier, or 
firn, allows to deploy antennas where temperatures are 
stable at— 52±2°C throughout the year and they are not 
affected by occasional high winds. A number of designs 
have been studied for this environment [5] using NEC4 
simulations. It was found that the fat tube-dipole, shown 
in Figure 1, offers a low power reflectivity in the fre- 
quency range above 25 MHz as well as a low group delay, 
matching the predictions from NEC4 simulations well 
(c.f. Figure 2). Manufactured from 5 cm copper tubes 
with a shrinking diameter towards the feedpoint, this de- 
sign has proven very robust for the six antennas of this 
type deployed at the South Pole thus far. 



2.2. Experimental setups 

To study the radio environment for air-shower detec- 
tion, two experimental setups have been installed in the 
2010/201 1 and 201 1/2012 polar seasons. Co-deployment 
with the ARA installations, and use of their DAQ system, 
allowed for a minimum investment. In the first season, 
two fat tube-dipole antennas have been deployed at a dis- 
tance of 30.5 m and at a relative azimuthal orientation of 
21 °. Both antennas were buried in trenches at a depth of 
~0.5m, and the trenches were backfilled with snow. The 
antennas are connected to the ARA-testbed DAQ sys- 
tem [6], where they are read out with a sampling rate of 
up to 2GS/s at 12 bit and a bandwidth of 30-300 MHz. 
In this setup, no dedicated trigger is available for these 
two surface antennas, so only events from a forced trig- 
ger every 2 s enter the analysis. In the second setup, four 
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FIGURE 2. Power reflectivity as a function of frequency for 
a fat tube-dipole from lab measurements and simulation. 



fat tube-dipole antennas were deployed, approximately 
in an equilateral triangle with one antenna in the cen- 
ter. Three of these antennas are aligned east- west, within 
±3 ° of the magnetic field, while one antenna is offset 
to that direction by +22 ° . To avoid different signal prop- 
agation times, all antennas are connected to the central 
DAQ by 200m of highly shielded cables and an overall 
gain of 74 dB. In this setup the signal from the four sur- 
face antennas in the 25 MHz- 120 MHz band is diplexed 
with the signal above 120 MHz from the V-pol antennas 
deployed in the ice for readout. In addition to the in-ice 
and minimum bias triggers, a dedicated surface trigger 
has been operated in a 3-out-of-4 mode with individual 
thresholds adjusted per channel. Due to low-temperature 
instability of the commercial single-board computer em- 
ployed in the DAQ system, only a limited amount of data 
has been collected with this setup from Jan, 26 th to May, 
3 rd 2012, with an uptime of less than 50%. 



2.3. Noise measurement 

Figure 3 shows the power spectral density, averaged 
over all waveforms, obtained in one hour as a function 
of the day from the forced trigger data in 2011. The 
signal is dominated by the falling spectrum of galac- 
tic noise from the low-pass boundary at 25 MHz up to 
about 90 MHz, above which it is limited by thermal noise 
in the DAQ chain. A single notch filter at 46.4 MHz is 
used to mitigate the signal of the meteor radar [7]. The 
absence of any other features and the very low over- 
all levels reflect the stability of the conditions at South 
Pole and the absence of any steady-mode radio transmit- 
ters. Figure 4 shows the power spectral density averaged 
over one hour for the surface triggered data obtained in 
2012. The enhancement of the galactic noise band and 
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FIGURE 3. Signal power spectrogram as a function of days 
in 2011 for minimum bias events taken with the ARA-testbed 
setup. 



the absence of contributions in the frequency range above 
50 MHz indicates that most triggers are caused by fluc- 
tuations in the thermal noise. Broadband emission not 
consistent with galactic noise is only observed for three 
short time ranges, possibly associated with exceptional 
air-craft communication activity. This absence of tran- 
sient broadband emitters makes South Pole a unique lo- 
cation for the radio detection of air- showers. 

We have also studied the correlation of the noise with 
wind speed, air pressure on the ground and air tempera- 
ture and found no significant dependence. A dependence 
of the trigger rate on the wind speed, as previously ob- 
served with radio antennas deployed with the RICE setup 
[8], has not been confirmed. This supports the hypothesis 
that these arise from discharge events on structures above 
the firn surface that are absent in the setup discussed here. 

To establish the overall calibration of the system we 
have used the galactic noise model by de Olivera [9], 
that provides interpolated full-sky intensity maps in the 
range from 10 MHz to 10 GHz. Convolving these maps 
for each frequency with the antenna response taken from 
the NEC4 simulation described in section 2.1, we obtain 
a noise spectrum prediction. Noise values for each com- 
ponent in the DAQ system are taken from manufacturer 
information and are traced through the full chain. Above 
25 MHz, the resultant spectrum agrees to within ±2dB 
with the measured spectrum. Since the antenna sensitiv- 
ity is not isotropic in azimuth, this agreements only holds 
if averaging the noise spectra over time periods much 
longer than one day. By convolving the galactic noise 
prediction with the antenna response instead for a partic- 
ular antenna orientation, i.e. a particular time of the day, a 
time dependant noise level is expected in each frequency 
band. Figure 5 shows the variation in noise RMS ampli- 
tudes for three frequency bands normalized to the daily 
average in each band. Higher noise levels are predicted 
in the lower frequency band when the antenna is aligned 
perpendicular to the direction of the galactic center (for 
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FIGURE 4. Signal power spectrogram as a function of days 
in 2012 for surface triggered events taken with the ARA-01 
setup. 



co-polarized radio- waves). However, for the highest fre- 
quency band, the largest noise levels are predicted when 
the antenna is oriented towards the galactic center (in 
cross-polarized mode). The inversion point, defined as 
the frequency at which the antenna becomes more sen- 
sitive in cross-polarized than in co-polarized mode, as 
well as the amplitude of the azimuthal modulation of 
the antenna response, function will depend on the elec- 
trical permittivity of the surrounding snow. From a X 2 
fit procedure we find the best agreement with the mea- 
sured noise level variations for a permittivity of £ = 1 .4 
at # 2 /NDF = 1.31. Different approximations have been 
proposed for the dependence of the permittivity of snow 
on the snow density [10, 11], suggesting a correspond- 
ing snow density of p S now = 

0.21 -0.29 g/cm 3 . Values of 
Psnow = 0.4g/cm 3 have been measured for pristine snow 
at the South Pole at a depth of one meter, but lower val- 
ues are plausible for the backfilled trenches in which the 
antennas are located. 
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FIGURE 5. Predicted (lines) and measured (dots) normal- 
ized noise RMS amplitude variation for three different fre- 
quency ranges as a function of the time of the sidereal day. 
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3. PERFORMANCE ESTIMATE 

Given the ideal environmental conditions outlined above, 
a Radio Air-Shower Test Array (RASTA) has been pro- 
posed to establish the feasibility of radio air-shower de- 
tection and study the physics potential [1]. A dedicated 
simulation has been developed to study the response. The 
proposed detector geometry, shown in figure 6, consists 
of 74 antennas arranged in pairs in a hexagonal grid with 
a spacing of 62.5 m, half the grid size of the IceTop array. 
In each pair, one antenna is oriented north-south and the 
other east- west relative to the geomagnetic field. 



3.1. Simulation chain 

Three different emission models have been used 
as inputs for the simulation chain: MGMR [12], 
REAS3.0 [13] and CoREAS [14]. For each model, we 
have simulated several hundred events in the energy 
range between 3 PeV to 1 EeV and with zenith angles up 
to 60°. The electric field as a function of time is sampled 
at each antenna pair location and folded with the antenna 
response (c.f. section 2.1). Galactic noise, following 
the model of Cane [15], and thermal noise are added 
assuming a feedpoint impedance of 50 £1. We simulate 
a second order bandpass filter from 25-300 MHz and 
assume digitization with 300 MHz at 16-bit resolution. 
An event is defined to trigger the array if 4 antennas in 
4 distinct pairs have a signal 5 times above the thermal 
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FIGURE 7. Trigger efficiency as a function of energy for 
proton (green) and iron (blue) primaries. 



noise within a time window of 1 /is. A sample trace 
of the simulated signal after digitization is shown in 
Figure 8. Figure 7 shows the trigger efficiency for proton 
and iron primaries from the CoREAS simulation which 
includes refractive index variation throughout the atmo- 
sphere. A 50% efficiency is reached at ~ 100 PeV, while 
for the REAS3.0 simulation, that assumes a constant 
refractive index, we find a 50% efficiency at ~ 30 PeV. 
It should be noted that, due to the strong dependence of 
the signal strength on the orientation w.r.t to the Earth's 
magnetic field and the zenith angle, a fraction of the 
events will also be detected at energies well below this 
threshold and with a steeply falling cosmic ray spectrum 
that may dominate the accumulated data. 
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FIGURE 6. Proposed layout of the Radio Air-Shower Test 
Array (RASTA) 



3.2. Event reconstruction 

In a first attempt of event reconstruction, we employ an 
iterative plane wave fit procedure. We perform a Hilbert 
transform to the signal in each antenna and fit a Gaussian 
to the peak (c.f. Figure 8). The peak time and amplitude 
of this Gaussian fit are used as a proxy to define a sig- 
nal hit if the amplitude exceeds five times the RMS am- 
plitude in the trace. Using singular value decomposition 
(which is faster and more accurate than x 2 minimiza- 
tion), we obtain the best plain wave best match to the set 
of antenna positions and arrival times. Excluding one of 
the single antenna hits at a time, we find the hit with the 
largest pull on the reconstructed direction and recursively 
remove these hits until the resulting direction does not 
vary by more than 0.3 ° (or not enough hits are left to con- 
tinue the procedure). Figure 9 shows the resulting angu- 
lar resolution for CoREAS simulated events. A median 
angular resolution of 1 .5 ° can be achieved even with this 
simple approach, comparable to the resolution of IceTop. 
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FIGURE 8. Simulated signal traces in one antenna pair from 
REAS3.0 including noise and system response. Time and am- 
plitude are extracted from a Gaussian fit (dashed line) to the 
Hilbert envelope. 
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FIGURE 9. Angular resolution from a plane wave fit that 
iteratively rejects hits with large pull factors. Iteration will stop 
if the fit converges (blue) or not enough hits remain (green). 



We also use the amplitude weighted geometric mean of 
all hits as a proxy for the shower core location, and find 
an average resolution of ~ 10m for REAS3.0 simula- 
tions not including refractive index variations. This sim- 
ple proxy does not work for CoREAS simulation. With 
refractive index effects included, the lateral distribution 
no longer falls monotonically with the radius for small 
distances to the shower axes [16]. 



tive signal cleaning algorithms, suggest that angular res- 
olutions of about 1 ° and shower core reconstructions of 
10 m may be obtained. This performance is comparable 
to the performance that has been measured for the Ice- 
Top detector. While successfully funded by the BMBF 
in 201 1, funding for logistic support for such an installa- 
tion at the South Pole has remained elusive, temporarily 
delaying future research for this promising project. 



4. SUMMARY AND OUTLOOK 

A series of experimental setups has been installed at the 
South Pole to explore the suitability of the location for 
an air- shower detector based on radio technology. Noise 
measurements from more than a year of data-taking 
demonstrated the excellent suitability of the South Pole 
for this approach. A detailed comparison of the measured 
noise spectra show agreement with model predictions on 
the few dB level, and hence a profound understanding of 
the antenna response and DAQ chain calibration. Moti- 
vated by this success, proposals have been submitted to 
the BMBF and the NSF for an extended radio air- shower 
test array — RASTA — consisting of 74 antennas de- 
ployed over a few 100 m 2 area. Performance estimates 
using particle-level shower simulations, as well as de- 
tailed simulation of the galactic noise, the DAQ response 
and antenna response show that a 50% detection effi- 
ciency at ~ lOOPeV may be achievable in such a geom- 
etry for a simple threshold-based trigger. Using straight- 
forward plane- wave reconstruction methods, with itera- 
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